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ABSTRACT 

We have obtained spectroscopic observations of the nebulae around seven Wolf- 
Rayet (WR) stars, four in the Milky Way (WR 8, 16, 18 and 40) and three in the 
LMC (BAT99-2, -11 and -38). They were observed using the ESO NTT with the 
EFOSC-2 instrument, except for one nebula, NGC 3199, which was observed using 
the UVES echelle spectrometer on the VLT. The aims of these observations were to 
(a) quantify the degree of chemical enrichment in WR nebulae which had previously 
been suggested to have abundances reflecting nucleosynthetic processing; and (b) to 
attempt to detect the far-red lines of neutral carbon (e.g. [C i] 9850 A) in the nebulae 
around WC stars. Nebular densities, temperatures and elemental abundances were 
derived using standard emission line diagnostics. For those nebulae for which the usual 
temperature diagnostic lines were not detected, their measured [N ii] and [O ii] line 
intensities were used to derive reliable N/0 ratios. Our high spectral resolution UVES 
datasct for NGC 3199 allowed the determination of line broadening temperatures 
using lines from several different species. These showed consistent patterns and fair 
agreement with the nebular temperatures derived from diagnostic line ratios. Amongst 
the Galactic WR nebulae, the nebulae around WR 8, 16 and 40 were found to be 
strongly nitrogen-enriched. NGC 3199 displayed abundances which were very similar 
to those of Galactic H ii regions, in agreement with previous analyses and consistent 
with the inference that it is a shell of swept-up interstellar gas. Amongst the LMC 
WR nebulae, none displayed an N/0 ratio that significantly exceeded LMC H ii region 
values. The far-red [C i] lines were detected from NGC 3199, with the strength of the 
8727A line indicating that at least part of it was produced by C"*" recombination 
rather than coUisional excitation. The line widths of the [C i] 9824 and 9850A lines 
were however the same as those of coUisionally excited nebular lines, indicating that 
those lines originated largely from ionized regions within the nebula. 

Key words: stars: Wolf-Rayet, stars: winds, outfiows, stars: circumstellar matter, 
stars: mass- loss 



1 INTRODUCTION 

Massive stars are known to be prodigious recyclers of the 
ISM, for example a 60 M© star is thought to end its life 
as a black ho le with a mass of around a few solar masses 
l|Smarttll2009l ). implying that that the star will expel around 
90% of its initial mass into the ISM during its lifetime. This 
mass loss, either as a smooth wind or a violent episode, 
can create circumstellar nebulae. We can help to quantify 
massive star recycling, in terms of its nucleosynthetic effects, 
by spectroscopically investigating the nebulae they create. 
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Spectroscopic observations of WR nebulae have been 
performed several times in the past, usually f ocusin g 
on the bright es t Gal ac tic ex a mples e. g. iKwitteJ (Il984h: 
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ex a mples e.g; 

1990l. Il99ll. Il992h : festeban fc Vilchej 
}l999i r~ As summarised by 
Barlow, (i201Ql . the nebular abundances found 



lEsteban et al 
(|l992l); 'Marston et alj 



tend to fall into three groups: no enrichments (wind blown 
bubbles, denoted W type nebulae); helium enrichment 
(wind blown bubbles with some degree of enrichment, 
denoted W/E type nebulae) and full blown He and N 
enrichments which are taken to be evidence for stellar 
ejecta (ejecta or E type nebulae). 

Massive stars have been mooted as possible sites of sig- 
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nificant carbon formation (|Henrvet al.ll2000l ). Wolf-Rayet 
stars in their WC phase possess very carbon rich atmo- 
spheres and high mass loss rates. Their role in enriching 
carbon has proven difficult to quantify given the challenge 
of observing gaseous carbon in nebulae around evolved stars. 
The coUisionally excited lines (CELs) of carbon are elusive, 
lying in relatively inconvenient parts of the electromagnetic 
spectrum: e.g. the [C ll] fine structure line in the far IR 
(158^m), various lines in the UV (e.g. C ill] 1909A) and 
[C l] hues in the far red (8700-9900A). Observations of the 
far-IR and UV lines rely upon satellite missions like Her- 
schel and the HST (Hubble Space Telescope). In contrast, 
the [C l] lines can be accessed from the ground, should they 
have sufficient intensity. 

The far-red carbon lines ([C l] 8727, 9824, 9850A) are 
routinely observed in certain astrophysical contex ts, partic- 
ularly from planetary nebulae and If ll regions. iLiu et ahl 
l| 19951 ) observed the lines in PN spectra, detecting all three 
lines in four nebulae. The lines are also common in the 
pho todissociated regions (PDR s) surrounding H ll regions, 
e.g.lMunch fc HippeleinI l|l982l l for NGC 2024. lBurton et al.1 
(| 19921 ) for NGC 2023. However it was subsequently shown 
that in these PDRs the lines must arise from recombina - 
tion rather than coUisional excitation ()Escalante et al."l99l|) 
which is also the case for some PNe fDanzigcr & Goad 1973) 
althoug h not for the PNe o bserved bv Xiu et al., (2995). 

Stock fc Barlowl l|2010l ) recently surveyed and classified 
nebulae around WR stars, discovering a new nebula around 
the WN/WC type star WR 8 and identifying an LMC WC 
star whose nebula would be appropriate for an abundance 
analysis. These nebulae, along with a "control group" of 
other WR nebulae were identified and observed with ESO's 
3.6m NT T/EFOSC2 (ESQ faint object spectrogra ph and 
camera) (|Buzzoni et al.lll984l : ISnodgrass et al.ll200^ ') in long 
slit spectroscopy mode with grisms covering regions from 
3100-lOOOOA that included a variety of standard nebular 
diagnostics. Complementary service mode observations of 
NGC 3199 around WR 18 were performed concurrently us- 
ing ESO's UVES instru ment mo unted on UT2 of the Very 
Large Telescope (VLT) (|Dekker ct al. 2000). 

From the data obtained we derive nebular physical con- 
ditions using standard nebular density and temperature in- 
dicators and then use these to derive elemental abundances 
from the CELs. We then compare these quantities with pre- 
viously published physical conditions and abundances for 
these nebulae and with those predicted by stellar evolution 
models. For the nebulae which are previously unobserved, we 
estimate whether there is any stellar ejecta content within 
the nebular material. 



2 TARGET SELECTION 

The LMC WC star BAT99-11 and the Milky Way WN/WC 
star WR 8 are carbon rich WR stars with ejecta-type cir- 
cumstellar nebulae. One of the aims of our observing pro- 
gramme was to attempt to detect the far-red carbon lines 
([C l] 8727, 9824, 9850A) in the nebulae surrounding both 
stars. Four other nebulae, each around a nitrogen rich (WN 
type) WR star, were chosen as a representative sample 
against which comparisons could be drawn; both with the 
nebulae with WC central stars and with nebulae analysed 



by previous authors. The parameters of each of the chosen 
targets are summarised in Table [1] 

Four of the targets (WR 16, NGC 3199, RCW58 and 
BAT99-2) have been the subjec t of previous spectrosc opic 
analyses, as briefly discussed bv lStock fc Barlow! l|2010l ). In 
general helium and nitrogen overabundances were found, 
however in each case the origin of the ejecta is still subject to 
debate in that it is possible that the nebular material could 
have been ejected during the red supergiant (RSG) or lu- 
minous blue variable (LBV) phases through which massive 
stars pass prior to becoming a WR star. 

The final target, the nebulosity thought to be associated 
with BAT99-38, has not previously been observed spectro- 
scopically. BAT99-38 resides in the LMC superbu bble N51D, 
and t herefore has quite corn plex surroundings (|Chu et al.l 
l2005h . IStock fc Barlow! (|2010l ') identified the arc of nebulos- 
ity near BAT99-38 as a possible ejecta nebula despite its 
presence within a superbubble. 



2.1 Previous Spectroscopic Observations 

RCW 58 was the subject of several spectroscop ic investiga - 
tions in the years following its identifi catio n bvlCh^ (|l98ll ) 
as a potential WR ejecta nebula. Kwitteil (|l984l ) presented 
the first long slit spectra of the nebulosity, which was found 
to be highly enriched in helium and nitrogen at an assumed 
temperature of 7500 K. Subsequ ently the dynamica l struc- 
ture of the nebula was revealed bv lSmith et al.l (| 19881 ) to be a 
combination of a quickly expanding shell interspersed with 
higher density clumps. It was speculated that the higher 
density material was the product of pre-WR stellar evolu- 
tion and did not represent ejecta from the WR phase. 

The anonymous circu mstellar nebula surrounding WR 
16 has been observed by iMarston et al. I l|l999t ) with sur- 
prising results. The nebular spectrum was devoid of several 
exp ected lines e.g. [O ill] 4959, 5007A; [S ii] 6717, 673lA. 
M arston et al.l (|l999l ) suggested that this was due to high 
densities (rie > 10'' cm~^ ) suppressing these lines. They 
noted that the [S ll] 6717, 673lA lines reappear in the outer 
regions of the nebula, which was assumed to indicate a much 
lower density of a tenth of the value which they assumed for 
the inner nebula. 

NGC 3199 has been the subject of several spectral 
analyses, albeit none with compara ble resolution to our 
VLT/UVES dataset. lEsteban etal] (1992) found abun- 
dances generally in accordance with those of Galactic H ll 
regions. 

The circumstellar neb ula around BAT99-2 wa s first ob- 
served spectr oscopically bv'Carne tt fc Chul (|l994h and sub- 
sequently bv lNaze et al . (2003) (the star was referred to in 
both cases as Brey-2) . Both found that the nebula was highly 
ionized - displaying He ll and [Ar iv] lines which are not 
norm ally observed in WR nebulae. BAT99-2 is a WN2 type 
star (|Foellmi et al.ll2003l ). amongst the hottest WR stars, 
with a surface temperature above 90,000 K (Crowthcr 200^), 
producing very intense EUV radiation and ionizing its sur- 
roundings to a greater degree than the nebulae around other 
WR stars. No significan t abundance enha ncements were de- 
tected in the nebula bv lNaze et al] l|2003l ). 
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Table 1. Target Parameters 



Nebula 


Central Star 


Spectral 
Type" 


V mag.'^ 




EfB-V)" 


RA(2000) 


Dec(2000) 


Galactic 
















Anon 


WR 8 (HD 62910) 


WN7/WC4 


10.48 


2.64 


0.85 


07 44 58.2 


-31 54 29.6 


Anon 


WR 16 (HD 86161) 


WN8I1 


8.44 


2.05 


0.66 


09 54 52.9 


-57 43 38.3 


NGC 3199 


WR 18 (HD 89358) 


WN4 


11.11 


2.92 


0.94 


10 17 02.3 


-57 54 46.9 


ROW 58 


WR 40 (HD 96548) 


WN8I1 


7.85 


1.56 


0.50 


11 06 17.2 


-65 30 35.2 


LMC 
















Anon 


BAT99-2 (Brey 2) 


WN2 


16.22 


< 0.5 




04 49 36.1 


-69 20 54.5 


Anon 


BAT99-11 (HD 32402) 


WC4 


13.95 


< 0.5 




04 57 24.1 


-68 23 57.3 


Anon 


BAT99-38 (HD 36402) 


WC4+0 


11.50 


< 0.5 




05 26 03.9 


-67 29 57.0 



a: from Ivan der Huchtl l|200lh (Galactic). [Brevsacher et al.l l|l999f l (LMC) 



3 OBSERVATIONS 
3.1 NTT/EFOSC2 

Most of the observations were performed in visitor mode 
over five days in December 2009 using the EFOSC2 instru- 
ment mounted on the 3.6m ESQ New Technology Telescope 
(NTT) at La Silla Observatory, Chile. Two grisms were used 
in long slit spectroscopy mode providing a wavelength cov- 
erage from 3095-5085A and 6015-10320A. The EFOSC2 slit 
length was 4.1 arcmin which, due to the fact that almost all 
of our targets had smaller angular sizes than this, enabled us 
to use the off-target portions of the slit for sky-subtraction. 
In all science observations the slit width was set to 1 arcsec- 
ond, while for standard stars the widest slit (5 arcseconds) 
was used. 

Grisms 14 (hereafter Blue) and 16 (hereafter Red) re- 
spectively have resolutions (FWHM) of SA and 14A respec- 
tively. This allowed us to use 2*2 binning throughout our 
observations without losing any information, as the line pro- 
files were still adequately sampled in this binning mode. This 
translated to 1030 pixels in both the dispersion and spatial 
axes. 

A log of the nebular observations is presented in Ta- 
ble [21 The slit positions used for each nebula are shown in 
Figures [THZ] A summary of the total integration times in 
both blue and red grisms is shown in Table [3] 

Multiple spectroscopic standard stars were observed 
each night, notably LTT 1788, HR718, and LTT 3864. Ob- 
serving templates for these stars existed in the EF0SC2 
database and were used each night so as to acquire the ap- 
propriate signal to noise level for each object to allow later 
fiux calibration. 



3.2 VLT/UVES 

Service mode VLT/UVES observations of NGC 3199 were 
peformed almost concurrently with the NTT/EF0SC2 ob- 
servations (see Table [2]). These were designed to assist with 
the main science goal while also potentially detecting heavy 
element recombination lines. The data cover the wavelength 
region from 3000-10400A with only small gaps between 
5700-5800A and 8500-8650A due to small gaps between 
CCD chips. The UVES dataset had a resolving power of 
^ — 40000 with a 0.6 arcsecond slit width and a 15 arcsec- 



Table 3. Total integration times by target 



Nebula 


Central 
Star 


Total Red 
(hours) 


Total Blue 
(hours) 


Anon 


WR 8 


4.5 


2.33 


Anon 


WR 16 


1.5 


1.5 


NGC 3199 


WR 18 


1.6 


1.6 


RCW 58 


WR 40 


1.5 


1.66 


Anon 


BAT99-2 


2 


2 


Anon 


BAT99-11 


6.33 


3.16 


Anon 


BAT99-38 


0.75 


0.75 



ond slit length. Its high resolution allowed us to forego the 
sky subtraction process, as night-sky emission lines appear 
at the instrumental resolution as opposed to the nebular 
lines for which we could measure their true widths. 

3.3 Individual Objects 

BAT99-2 

The slit position adopted for the LMC WR nebula around 
BAT99-2 is shown in Figure [1] This position was adopted as 
it placed the greatest extent of the arc shaped nebulosity un- 
der the EF0SC2 slit. The slit position we adop ted roughly 
matches the position described as "arc" by iNaze et al.l 
(2003). 

BAT99-11 

The slit position adopted for the nebula around BAT99-11 
in the LMC is shown in Figure [2] The nebulosity around 
BAT99-11 is physically much larger than that normally as- 
sociated with WR ejecta nebulae, so the slit position was 
chosen to intersect the inner portions of the nebula as this 
may be more recently expelled material. 

BAT99-38 

The slit position adopted for the nebulosity near BAT99-38 
is shown in Figure [3] It was discovered later that BAT99- 
38 was not the star which appears at the focus of the arc, 
but is embedded in the arc itself. This was due to a small 
error in the WCS system of the digitised SHS survey plates 
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Table 2. Journal of Nebular Observations 



Nebula 


Central Star 




Date 


Instrument 


Grism 


AA (A) 


Exposure 
Time (s) 


No. Ex- 
posures 


Anon 


vv rx, o 


2009 Dec. 


19/20 




14 




1 son 


2 


Anon 


vv IX o 


2009 Dec. 


19/20 






UU-LO--LUOZU 


1 son 




Anon 


WR S 
vv rx o 


2009 Dec. 


20/21 




14 




1 900 




Anon 


vv rt o 


2009 Dec. 


20/21 




ID 


UUIO-IUOZU 


1 f^OO 


D 


Anon 


WR 1 f\ 

vv IX 1\) 


2009 Dec. 


21/22 




14 




1 soo 


Q 
O 


Anon 


vv -tx ID 


2009 Dec. 


21/22 




ID 


uUlO-lUtjZU 


1 HOO 


O 




WR 1 S 

vv -TV ±o 


2009 Dec. 


28/29 


TT^/F'^ 

U V H/O 


DICl 


OUOU-OOOU 4 / UU-DO^U 


1 4Q0 




l\TPn "il QQ 


vv rt lo 


2009 Dec. 


28/29 


TT\/T?Q 




o / ou-^yyu dduu-iuduu 


1 zlQO 


9 


NPP "^l QQ 


WR 1 S 

vv -TV, lo 


2010 Jan. 1/2 


U V H/o 


DIC2 


4QQn fifion i of^no 
o ( ou-^yyu o6 uDuu-ruuuu 


1 4Q0 


2 


i\w vv oo 


WR An 

vv IX 4U 


2009 


Dec. 


18/19 




14 


0V170-0V00 


fion 

uuu 




RCW 58 


WR 40 


2009 


Dec. 


18/19 


EFOSC2 


16 


6015-10320 


600 


5 


RCW 58 


WR 40 


2009 


Dec. 


19/20 


EFOSC2 


14 


3095-5085 


600 


5 


RCW 58 


WR 40 


2009 


Dec. 


19/20 


EFOSC2 


16 


6015-10320 


600 


5 


Anon 


BAT99-2 


2009 


Dec. 


22/23 


EFOSC2 


14 


3095-5085 


1200 


3 


Anon 


BAT99-2 


2009 


Dec. 


22/23 


EFOSC2 


16 


6015-10320 


1200 


3 


Anon 


BAT99-2 


2009 


Dec. 


21/22 


EFOSC2 


14 


3095-5085 


1200 


3 


Anon 


BAT99-2 


2009 


Dec. 


21/22 


EFOSC2 


16 


6015-10320 


1200 


3 


Anon 


BAT99-11 


2009 


Dec. 


18/19 


EFOSC2 


14 


3095-5085 


1500 


4 


Anon 


BAT99-11 


2009 


Dec. 


18/19 


EFOSC2 


16 


6015-10320 


1500 


8 


Anon 


BAT99-11 


2009 


Dec. 


22/23 


EFOSC2 


14 


3095-5085 


1800 


3 


Anon 


BAT99-11 


2009 


Dec. 


22/23 


EFOSC2 


16 


6015-10320 


1800 


6 


Anon 


BAT99-38 


2009 


Dec. 


19/20 


EFOSC2 


14 


3095-5085 


900 


3 


Anon 


BAT99-38 


2009 


Dec. 


19/20 


EFOSC2 


16 


6015-10320 


900 


3 




Figure 1. SHS Ha image l lParker et al.ll2005l) showing the loca- 
tion of the slit for the nebula around BAT99-2. Entire slit length 
(5') is not shown. 




« i 



Figure 2. SHS Ha image showing the location of the slit for the 
nebula around BAT99-11. Entire slit length (5') is not shown. 



into question our preliminary morphological identification of 
this nebula as a possible ejecta nebula. 



which were used to id entify WR nebulae in the LMC (see 
IStock fc Barlowl (|201(]| )). The slit position therefore unex- 
pectedly included BAT99-38 along with the nebulosity. The 
slit was subseqently moved slightly such that BAT99-38 was 
no longer included in the slit. The position of the star called 



WR 8 

As with BAT99-11, it is unclear whether the knots which 
appear interior to the outer ring of nebulosity are in the 
same shell as the visible ring, or are really spatially closer to 
WR 8. The slit position was chosen to encompass the two 




Figure 3. SHS Ha image showing the location of the slit for the 
nebulosity near BAT99-38. Entire slit length (5') is not shown. 
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Figure 4. SHS Ha image showing the location of the slit for the 
nebula around WR 8. Entire slit length (5') is not shown. 

brightest Ha emitting knots which did not coincide with 
stars. The slit position adopted is shown in Figure |4] 

WR 16 

The adopted sht position for WR 16 is shown in Figure [5] 
Due to the large size of this nebula the slit shown in Fig- 
ure[5]is appoximately the true length of the EFOSC2 slit (~ 
5'). The slit position was chosen such that it encompassed 
some of the brighter Hq knots seen in the SHS imagery, 
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Figure 5. SHS Ha image showing the location of the slit for the 
nebula around WR 16. 




Figure 6. SHS short red image of NGC 3199 (location of WR 18 
indicated with dashes), ^ 15 arcsecond VLT/UVES slit is shown 
in the expanded section. 



along with enough off-nebula sky to perform a reasonable 
sky subtraction. 



NGC 3199 (WR 18) 

The slit position adopted for NGC 3199 is shown in Figure|6l 
Given the shorter (15") UVES slit length the brightest knot 
visible in the short red SHS exposure was chosen. For this 
object the Ha SHS image is saturated and so not useful 
in terms of selecting the brightest parts of the nebulosity. 
While NGC 3199 is very much brighter in Ha than our other 
targets, it was still important to select the brightest possible 
section of nebulosity such that the likelihood of detecting 
weak lines was maximised. 
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Figure 7. SHS Ha image showing the location of the sht for 
RCW 58 (WR 40). 



RCW 58 (WR 40) 

The slit position used for observations of RCW 58 is shown 
in Figure [7] This sht position was adopted to maximise the 
amount of bright nebulosity on the western side of RCW 58 
in the slit. 



4 DATA REDUCTION 
4.1 EFOSC2 Spectra 

At the beginning and end of each night a series of calibra- 
tion frames were acquired, including bias, arc and flat field 
frames. The reduction of the data was performed using ESO- 
MIDAS l|Warmel j [1993 ) and Gaussian line measurements 
were performed using the ELF emission line fitting routine 
written by P.J. Storey within the Dipso Starlink package 
written by I. Howarth. A typical result of the data reduc- 
tion process, a 2D spectral image showing nebular lines, is 
shown in Figure |8] 

Sky subtraction was performed using the ESO-MIDAS 
routine "SKYFIT/LONG", which, given two sky windows 
either side of the object spectrum, calculates a 2D frame 
of the night sky lines by interpolating across between the 
windows. This can then be subtracted from the data frame 
to leave just the object spectrum in the resulting frame. 
The 7000-9000A region is particularly blighted by night sky 
lines, as shown in the composite Figure |5] which displays a 
non-sky subtracted frame, an interpolated sky frame and the 
sky-subtracted image. Post sky-subtraction, some evidence 
of fringing appears which we were not able to remove using 
any combination of internal flats or dome flats. This residual 
fringing is not thought to be a problem as there were very 
few lines being measured in the 7000-9000A spectral region. 

For each 2D co-added sky-subtracted nebular spectrum, 
every row of pixels which was thought to contain an un- 



Figure 8. 2D Spectrogram of the nebulosity around WR 8 in the 
3000-5050A wavelength region (Grism #14). The nebular lines at 
the top are from the nebulosity shown in Figure |4] The emission 
line source in the lower portion of this spectrum is the scattered 
light of the WR central star. 



contaminated nebular spectrum was co-added to produce 
ID spectra. Subsequently each line was measured using 
Dipso/ELF. 

The reddening was derived using the Balmer series ra- 
tios to calculate c(H/3) - the logarithmic extinction at H/3 - 
using pairs of Balmer lines (one of which was always H/3). 
The derived reddening values (Table[4|) were consistent with 
literatur e values derived by previous au thors in the ca ses of 
WR 16 (iMarston et al.lll999 '). WR 40 (|Kwitteil [l983 ) and 
BAT99-2 fNaze et al.ll2003l ). T hey were also checked against 
the Galactic extinction maps of lSchlegel et al.l(|l998l ). In ev- 
ery case the Galactic WR nebulae had derived E(B-V) val- 
ues lower than the Schlegel et al maximum value (indicating 
that they are, in fact, resident in the Galaxy - a useful san- 
ity check). However for the LMC nebulae the Schlegel dust 
maps indicated E(B-V) values which were too high by at 
least a factor of f i ve. Fo r LMC nebula we resorted to the 
iBurstein fc Heile j (|1982| ') E(B-V) maps which were found to 
agree much better with our data. 

In some cases the dereddening process highlighted a lack 
of a reliable relative flux calibration between the blue and 
red spectra. For the nebulae around WR 16 and BAT99-11, 
this manifested itself as very different c(H/3)'s when com- 
paring Ha in the red spectra to H/3 in the blue, compared 
with those derived from just the blue spectra (using H(5 and 
H7). For the nebula around WR 8, the converse situation 
applied in that a very low c(H/3) was found from the blue 
Balmer lines and a high c(H/3) was found from Ha. This may 
have been caused by changing conditions throughout our al- 
located nights, making the flux calibrations derived for the 
beginning and end of each night inappropriate for the mid- 
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Figure 9. left: 2D spectra of the region around WR 8 for tlie red spectral region (6000-lOOOOA Grism #16); centre: ESO-MIDAS 
"SKYFIT/LONG" generated 2D sky spectrum for the frame shown; right: Sky subtracted spectrogram of the region around WR 8. 



die of each night. This effect was minimised by using the 
nearest (in time) flux calibration to each observation. The 
adopted c(ff/3) for both blue and red settings, where they 
did not agree, are shown in Table [l] In cases where a sepa- 
rate reddening solution was adopted for the red specta, the 
calibration scales the red spectra such that Ho? has the cor- 
rect flux ratio to the Balmer lines in the blue spectra. This 
has the consequence of ensuring that the region around Ha 
also posesses the correct relative flux levels (e.g. [N ll], [S ll]). 
Spectral lines detected far from Ifa were considered unreli- 
able in these cases. 

A futher complication was intoduced by the choice to 
deviate from the parallactic angle in order to maximise the 
amount of nebulosity falling under the slit. In cases where 
the nebulae were large and of fairly constant surface bright- 
ness, this did not cause a problem. However, the nebulae 
with discordant blue/red reddening values were those where 
the slit passed through several small clumps, e.g. WR 8 (See 
Figure [4]). In these cases the effects of atmospheric disper- 
sion combined with the small angular scales of the clumps 
of nebulosity may be the cause of the reddening issues. The 
offsets from the parallactic angle for each object are listed 
in Table [1 

By adopting a different reddening solution for the red 
spectra we effectively de-redden and flux calibrate Ha and 
the adjacent lines ([N ll], [S ll], He I 6678A) but lose flux cal- 
ibration for lines far from Ha ([S ill] 9069, 953lA for exam- 
ple). For the nebulae which required this technique (those 
around WR 8, WR 16 and BAT99-11), this would render 
their [S in] diagnostic temperatures unreliable. 

Due to the relatively low resolution in both the blue and 
red spectra, blended lines could be a problem in these data. 
Most seriously, the standard density diagnostic line pairs 
[S ll] 6717, 673lA and [O ll] 3727. 3729A were blended to 
some degree. The [O ll] doublet is completely unresolved: its 
measured FWHM is very close to the 8A instrumental reso- 
lution. The [S ll] doublet is marginally resolved by the 14A 
resolution in the red spectra, in that there are clearly two 
components. This situation is far from ideal as this leaves 
the [S ll] doublet as the only density diagnostic, and with 



a larger uncertainty on the [S ll] doublet ratio due to its 
semi-blended nature. 

It was also found that the observed intensity ratio of the 
[O III] hues of 5007A to 4959A was not consistent with the- 
ory. This intensity ratio is fixed by the rela tive atomic transi- 
tion p robabilities and has a value of 2.98 (|Storev fc ZeippenI 
l2000t ). whereas we observed it to be nearer 2.0 in all cases. 
This was attributed to the 5007A line falling very close 
(within 3-4 pixels) to the edge of the detector, and the flux 
calibration being unreliable there. It was assumed that the 
4959A line was correctly flux calibrated and so it was used 
on its own at the [O ill] temperature diagnostic stage. 



4.2 UVES Spectra 

The calibrations for UVES observations that are provided 
by ESO are intended for point source objects, however, and 
included a sky-subtraction step which was inappropriate for 
our observations. By default the UVES reduction pipeline 
assumes that two strips along the edges of the slit are free 
of the science target spectrum and treats them as sky pix- 
els. The pipeline then averages between these sky pixels and 
subtracts this from the overall reduced spectrum. The neb- 
ulosity of NGC 3199 is very much larger than the ~ 15" 
slit and as such no region of the slit represented a "sky" 
spectrum. The pipeline was therefore re-run with the sky 
subtraction phase excluded for all the UVES observations. 

Cosmic ray subtraction was performed for the UVES 
observations by median combination of our four science 
frames for each setting. This produced a final science frame 
for each setting. The flux calibration was checked by com- 
paring line fluxes in the overlap regions between detectors, 
which were found to agree to better than 10%. 

The EF0SC2 problem of line blends was, of course, ab- 
sent from the UVES observations and lines which are usu- 
ally blended in lower resolution spectra are well resolved, 
e.g. [O ll] 7319, 7320, 7330, 733lA. In fact the higher reso- 
lution of the UVES instrument made the velocity structure 
of the nebula evident in the bright lines. The best resolu- 
tion of this velocity structure is seen in the far red [S ill] 
lines, due to their high atomic weight minimising thermal 
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Table 4. Derived Reddenings 



Nebula 



Central 



c{H/3) c{H/3)-^^ <^l^)red I-^Plied Stellar 



Maximum Literature 





Star 






E{B-V)' 


E(B-V)'= 


Galactic 
E(B-V)'' 


E(B-V) 


Galactic 
















Anon 


WR 8 




0.04 1.73 


0.02 


0.85 


1.23 




Anon 


WR 16 




0.99 2.43 


0.68 


0.66 


2.97 


0.44 


NGC 3199 


WR 18 


1.18 




0.81 


0.94 




0.79 
















1.10 


RCW 58 


WR 40 


1.07 




0.74 


0.50 


0.79 


0.77 


LMC 
















Anon 


BAT99-2 


0.42 




0.29 




0.09 


0.20 


Anon 


BAT99-11 




0.16 0.62 


0.11 




0.06 




Anon 


BAT99-38 


0.06 




0.04 




0.05 





Literature 
Ref. 



E(B-V) 



Marstonetal 
KwitteJ \l~A ) 



Estebanetal 
KwitteJ lll98l) 



lll999l ') 



lll992h 



iNaze et al.l ll2003h 



a: In some cases different c(H/3) w ere adopted for blue/red settings, see text. 
h: Calculated u sing relations f r om iHowartbl lll983l) with c(H/3) or c(H/3)5,„e 
c: Ivan der Huchd | |2001|'| Galactic; unknown LMC 
d: ISchlegel et al.l lll998l'l Galactic: iBurstein fc Heilej lll982l ') LMC 



broadening effects (Figure I10|) . It was found that the best 
fit to this structure utilised four components, two of which 
comprised the majority of the emission while a further two 
are distinctly seen redward of the main emission. 

The lines generated by lower mass ions are subject to 
a greater degree of thermal broadening, widening the red- 
shifted components to the point where they are indistin- 
guishable. As such it proved difficult to derive a separate 
reddening solution for the second velocity component. In 
fact it proved difficult to fit the line profiles with a consis- 
tent number of velocity components. The adopted solution 
was to fit the brighter component with one gaussian, if it 
produced an accurate fit (the strongest lines: Ha, [N ll], 
[O in] etc) and to use an extra component where it was 
required and subsequently to sum the components in the 
bright section. 

For NGC 3199 the reddening was derived from the 
UVES spectrum Balmer decrement. The c(H/3) that we de- 
rived, listed in Table [4] with those of the other n e bulae , 
largely agree d with prev ious derivations bv iKwitteJ (| 19841 ^ 
and Estcban et al.l 11993). The minor differences in redden- 
ing are likely due to the fact that each study has observed 
a different section of a large nebula. 



5 RESULTS 

Observed and dereddened relative line intensity lists for each 
nebula observed with EF0SC2 can be found in Tables [S] 
(LMC) and [6] (Galactic). A corresponding line list for NGG 
3199 is shown in Table [7] 

The [C l] 8727, 9824, 9850A lines were not de- 
tected in our NTT/EFOSC2 observations. However, in the 
VLT/UVES observations of NGC 3199 we did detect the 
[C l] lines, giving an example of what their expected inten- 
sity might be. The ratio of the dereddened intensities of [C l] 
9850A to [S III] 953lA in NGC 3199 is 0.006. If we apply 
this ratio to the NTT/EFOSC2 observations we find ex- 
pected line fiuxes of less than half the faintest detected line. 
We conclude from this that the [C l] lines are not strong 




Velocity, km s ' 

Figure 10. UVES velocity profile of [S in] 953lA for NGC 3199, 
observed profile in black, four-component Dipso/ELF fit dotted. 

enough in any of the NTT/EF0SC2 sample to be detected 
with the configuration which we utilised, or indeed with any 
of the available grisms for EF0SC2. 

The detection of [C l] 8727, 9824, 9850A in NGC 3199 is 
largely due to the high resolution of the UVES instrument. 
The enhanced resolution did not yield the detection of any 
heavy element recombination lines. The brightest carbon 
recombination line (C ll 4267A) was not detected despite 
experimenting with spectral binning of the high resolution 
observations to attempt to increase the S/N. The relatively 
high reddening towards NGC 3199 (E(B-V) = 0.81), does 
not help the detection of this blue line. 



6 ANALYSIS 

6.1 Plasma Diagnostics 

In total, four different density diagnostic line doublets were 
employed during the analysis. The [S ll] doublet was present 
in most cases, however the [O ll], [Ar iv] and [CI ill] doublets 
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Table 5. EFOSC2 Line Intensities (H/3 = 100) - LMC Nebulae 



Line 


Ion 


Observed 


± 


Dereddened 


± 


Observed 


± 


Dereddened 


± 


Observed 


± 


Dereddened 


± 








BAT99-2 






BAT99-11 






BAT99-38 




3727.0 


[O n] 


58.53 


0.72 


75.00 


0.92 


260.95 


0.97 


286.59 


1.06 


168.10 


8.94 


174.01 


9.26 


3868.7 


[Ne III] 


53.11 


0.69 


66.36 


0.86 


4.21 


0.80 


4.58 


0.87 


9.58 


6.91 


9.89 


7.13 


4101.7 


HS 


24.49 


1.43 


29.19 


1.71 


23.23 


0.80 


24.82 


0.85 


24.51 


0.70 


25. 12 


0.72 


4340.5 


H7 


42.11 


1.43 


47.60 


1.61 


46.15 


0.79 


48.33 


0.83 


43.66 


0.69 


44.41 


0.71 


4363.2 


[O m] 


12.85 


1.42 


14.45 


1.60 


2.76 


0.79 


2.88 


0.82 










4471.5 


He I 


1.77 


0.32 


1.94 


0.35 


4.73 


0.78 


4.90 


0.81 


3.82 


0.70 


3.87 


0.70 


4685.7 


He II 


72.52 


0.68 


75.61 


0.71 


















4711.4 


[Ar iv] 


13.21 


0.63 


13.69 


0.65 


















4740.2 


[Ar iv] 


8.64 


0.62 


8.89 


0.63 


















4861.3 


H/3 


100.00 


0.82 


100.00 


0.82 


100.00 


0.81 


100.00 


0.81 


100.00 


1.13 


100.00 


1.13 


4958.9 


[O III] 


212.87 


1.48 


208.00 


1.45 


140.47 


0.83 


139.25 


0.82 


76.92 


0.74 


76.67 


0.74 


5006.8 


[O III] 


592.93 


1.73 


572.82 


1.67 


352.67 


0.96 


348.10 


0.95 


178.59 


0.85 


177.73 


0.85 


6312.1 


[S III] 


2.62 


0.33 


1.99 


0.25 


3.06 


1.38 


2.04 


0.92 


0.71 


0.27 


0.69 


0.26 


6562.8 


Hq 


394.99 


1.64 


289.88 


1.20 


449.79 


3.05 


284.95 


1.93 


290.15 


2.07 


277.89 


1.98 


6583.5 


IN II] 


6.55 


1.42 


4.79 


1.04 


20.82 


2.51 


13.13 


1.58 


13.80 


1.69 


13.22 


1.62 


6678.2 


He I 










6.29 


2.51 


3.89 


1.55 


3.44 


1.71 


3.29 


1.63 


6716.4 


[S II] 


18.95 


1.66 


13.62 


1.19 


18.36 


3.15 


11.28 


1.93 


15.83 


2.12 


15.11 


2.03 


6730.8 


[S n] 


11.29 


1.67 


8.10 


1.20 


12.93 


3.12 


7.92 


1.91 


11.66 


2.12 


11.13 


2.02 


7005.4 


[Ar v] 


3.43 


1.30 


2.37 


0.90 


















7135.8 


[Ar III] 


11.16 


1.39 


7.61 


0.95 


19.45 


2.51 


11.05 


1.42 


9.72 


1.67 


9.21 


1.58 


9068.6 


[S III] 


16.10 


0.77 


9.14 


0.43 


29.56 


0.71 


12.84 


0.31 


10.98 


0.77 


10.15 


0.71 


9530.6 


[S III] 


22.26 


1.33 


12.25 


0.73 


49.91 


0.88 


20.69 


0.36 


5.03 


1.15 


4.62 


1.06 



Table 6. EFOSC2 Line Intensities (H^ = 100) - Galactic Nebulae 



Line 


Ion 


Observed 


± 


Dereddened 


± 


Observed 


± 


Dereddened 


± 


Observed 


± 


Dereddened 


± 








WR 8 






WR 16 






WR 40 




3727.0 


[O II] 


60.72 


6.98 


62.05 


7.13 










22.43 


6.32 


42.23 


11.91 


4101.7 


H5 


23.93 


4.04 


24.30 


4.11 


22.63 


2.12 


34.28 


3.22 


17.36 


0.86 


27.17 


1.35 


4340.5 


H7 


48.52 


4.03 


49.05 


4.07 


28.86 


2.12 


38.56 


2.84 


35.20 


0.87 


48.12 


1.19 


4471.5 


He I 














6.00 


0.85 


7.60 


1.08 






4861.3 


H/3 


100.00 


4.20 


100.00 


4.20 


100.00 


2.52 


100.00 


2.52 


100.00 


1.03 


100.00 


1.03 


4958.9 


[O III] 


109.15 


4.27 


108.93 


4.26 


7.17 


2.18 


6.79 


2.06 










5006.8 


[O III] 


226.13 


4.84 


225.45 


4.82 


7.67 


2.12 


7.07 


1.95 


8.82 


0.86 


8.08 


0.79 


6548.1 


[Nil] 


108.59 


7.83 


30.40 


2.19 


288.38 


14.40 


48.57 


2.42 


144.02 


5.77 


65.73 


2.63 


6562.8 


Ha 


1026.07 


9.40 


284.87 


2.61 


1711.12 


16.04 


284.82 


2.67 


631.23 


6.23 


286.59 


2.83 


6583.5 


[Nil] 


438.50 


4.91 


120.32 


1.34 


1140.75 


12.77 


186.78 


2.09 


543.30 


5.04 


244.89 


2.27 


6678.2 


He I 


43.01 


4.69 


11.19 


1.22 


55.69 


12.10 


8.46 


1.84 


17.91 


4.67 


7.81 


2.03 


6716.4 


[S II] 


16.99 


6.30 


4.32 


1.60 










13.82 


5.51 


5.95 


2.37 


6730.8 


[S II] 


11.27 


6.13 


2.84 


1.55 










9.92 


5.47 


4.25 


2.34 


7135.8 


[Ar III] 


47.69 


4.63 


9.76 


0.94 


43.66 


10.92 


4.74 


1.18 


13.58 


4.75 


5.11 


1.78 


9068.6 


[S III] 


65.30 


6.31 


6.28 


0.60 










19.45 


1.71 


4.60 


0.40 


9530.6 


[S III] 


84.69 


9.40 


7.14 


0.79 










41.63 


2.56 


9.07 


0.55 



were also used when tlie instrumental resolution and nebular 
ionisation permitted. 

All but one of the nebulae (WR 16) displayed the [S 11] 
doublet at 6716 A, 6731 A, which allowed calculation of their 
electron densities. Possible reasons for the absence of these 
lines in the nebula around WR 16 are discussed in Section[7l 
In each case where the [S 11] doublet was detected, its ratio 
indicated electron densities near the low density limit (1 < 
He <~ 300 cm~^). In cases where the [S 11] doublet ratio 
indicated a density less than 50 cm~"^, a density of 50 cm~"^ 
was adopted. As discussed in Section g] the NTT/EFOSC 
detections of the [S 11] doublet are blended. This results in 
a relatively high uncertainty on both the line measurements 
and subsequently the ratio of the two lines. In contrast the 
ratio of the [S 11] lines detected in NGC 3199 with UVES 
carries a much lower uncertainty as both lines were well 
resolved. 

The higher excitation nebulosity around BAT99-2 dis- 



played the [Ar iv] 47IIA, 4740A doublet. It yielded a density 
consistent with the [S 11] density, i.e. the low density limit. 

The increased resolution and wavelength coverage of the 
UVES instrument allowed the use of the other two line pairs 
mentioned, [O 11] and [CI in]. The densities yielded were con- 
sistent with each other and with the [S 11] density (see Ta- 
ble The [O 11] density in particular agrees very well with 
the [S 11] density and carries a low uncertainty. The [CI ill] 
lines are comparatively weak, so the uncertainty in the ratio 
is much higher. This results in a less well constrained diag- 
nostic, albeit one that peaks in roughly the same place as 
the [O 11] and [S 11] densities. 

Since the derived electron densities were all below the 
critical densities of the abundance diagnostic lines, the un- 
certainties associated with the electron densities have no 
effect on the abundance uncertainties we will derive. 

In general, the electron temperature diagnostics are 
more important to the final abundance determinations due 



10 D.J. Stock et al. 



Table 8. Nebular Diagnostics - NTT/EFOSC Observations 



Nebula 


[S 11] 

6716A 
6731A 


ne([S II]) cm-» 


[Ar iv] ne([Ar iv]) cm^-* 
4711A 
4740A 


[^"'] 
9069A+9531A " 

6312A 


Te([S III]) K 


[O III] 

4959A 
4363A 


Te([0 III]) K 


BAT99-2 

BAT99-11 

BAT99-38 


1.66 ± 0.29 
1.41 ± 0.43 
1.35 ± 0.31 


50'' 

t.„+440 

78t2f 


0.65 ± 0.06 50'' 


15 ± 3 

17tf 
41+f 


20000^= 
20000^= 
11200 ± 2200 


14.4 

48 


16300 ± 1000 
10400 ± 1000 



[s III] 9531A was tellurically absorbed, its intensity was calculated using the ratio I(953lA)/I(9069A) = 2.48 

'': Low density limit. 
High temperature limit. 



Table 7. UVES Line Intensities (H^ = 100) - NGC 3199 



Line 


Ion 


OK A 


-(- 


Dereddenedib 


3726.0 


[O n] 


87.62 


0.43 


175.51 


0.86 


3728.8 


[O II] 


107.42 


0.52 


214.88 


1.05 


3868.7 


[Ne III] 


36.24 


0.22 


67.62 


0.41 


3967.5 


[Ne III] 


11.16 


0.19 


19.72 


0.34 


4101.7 


m 


15.11 


0.22 


24.70 


0.36 


4340.5 


H7 


32.91 


0.25 


46.38 


0.35 


4363.2 


[O III] 


2.45 


0.08 


3.41 


0.11 


4471.5 


He I 


3.37 


0.13 


4.37 


0.17 


4861.3 


H/3 


100.00 


0.24 


100.00 


0.24 


4958.9 


[O III] 


218.60 


0.81 


204.88 


0.76 


5006.8 


[O III] 


680.69 


2.26 


618.02 


2.05 


5517.7 


[CI III] 


1.85 


0.14 


1.22 


0.09 


5537.6 


[CI III] 


1.40 


0.08 


0.91 


0.05 


5577.3 


[O I] 


0.25* 


0.22 


0.16 


0.14 


5875.7 


He I 


21.93 


0.14 


12.27 


0.08 


6300.3 


[O I] 


17.17 


0.39 


8.01 


0.18 


6312.1 


[S III] 


5.92 


0.07 


2.75 


0.03 


6363.8 


[O I] 


5.78 


0.14 


2.62 


0.06 


6548.1 


[N II] 


91.97 


0.40 


38.90 


0.17 


6562.8 


Ha 


675.20 


1.33 


283.96 


0.56 


6583.5 


[N II] 


287.83 


1.16 


120.09 


0.48 


6678.2 


He I 


8.79 


0.11 


3.54 


0.04 


6716.4 


[S II] 


118.35 


0.88 


46.96 


0.35 


6730.8 


[S II] 


96.39 


0.65 


38.04 


0.26 


7135.8 


[Ar III] 


61.93 


0.37 


21.19 


0.12 


7318.9 


[O II] 


3.06 


0.07 


0.99 


0.02 


7320.0 


[O II] 


11.40 


0.08 


3.68 


0.02 


7329.7 


[O II] 


6.07 


0.07 


1.95 


0.02 


7330.7 


[O II] 


6.15 


0.07 


1.98 


0.02 


8727.1 


[C I] 


1.00 


0.07 


0.22 


0.01 


9068.6 


[S III] 


480.20 


5.56 


98.72 


1.14 


9530.6 


[S III] 


1416.72 


14.39 


266.44 


2.70 


9824.1 


[C I] 


3.44 


0.20 


0.62 


0.03 


9850.2 


[C I] 


8.94 


0.13 


1.59 


0.02 



* upper limit 



to the Boltzmann factor dependence of the emissivity of col- 
hsionally excited lines. With this in mind, it was of impor- 
tance to exploit all the available temperature diagnostics. 
The NTT/EFOSC observations were designed to access the 
[O III] 5007, 4959, 4363 A Unes in the blue and the [S in] 
9532, 9069, 6312 A lines in the red. These diagnostics are 
strongly dependent upon the very weak, temperature sen- 
sitive, [O III] 4363 A and [S lll] 6312 A hues. These lines 
were detected in some cases, typically where high tempera- 
tures might be expected (BAT99-2 for example). The [S ill] 
6312 A line detections, in particular, should be treated as 
upper limits as they were typically blended with [O l] 6300. 



"D 
0) 




-10 10 
Velocity, km s' 

Figure 11. Comparison of the NGC 3199 flux normalised veloc- 
ity profiles of [C l] 9850A (dot dash), 9824A (dashed) and 8727A 
(solid). The [C l] 8727A line appears to have a narrow recombi- 
nation component. 



The VLT/UVES observations yielded a richer set of 
temperature diagnostics. In addition to the [O ill] and [S ill] 
lines, which should be regarded as the primary temperature 
diagnostics, temperature diagnostics involving [O l], [O 11] 
and [C l] were also detected. The UVES instrument, in the 
dichroic configuration we utilized, has an unfortunate wave- 
length coverage gap between the red and blue arms of the 
spectrograph in the region of [N 11] 5755 A which precluded 
the use of the [N 11] 6584, 6548, 5755 A lines as a tempera- 
ture diagnostic. The [O 11] temperature diagnostic (7319A 4- 
7320A + 7330A + 7331A versus 3726A + 3729A) was calcu- 
lated but was ultimately disregarded for two reasons. Firstly, 
the uncertainty on the diagnostic ratio is rather high due to 
the weakness of the [O 11] lines at 7325A; secondly, this ratio 
has proven vulnerable to several further systematic uncer- 
tainties due to the large wavelength separation between the 
lines and the possibility of a significant re combination lin e 
contribution to the 7320-1-7330 A lines (e.g. lLiu et"aLll200Gl '). 



6.2 The [C l] far-red lines 

The velocity profiles for the [C l] lines as observed in NGC 
3199 are shown in Figure[TTl The [C l] 8727A line appears to 
have two distinct components, one narrow and one of similar 
width to the other [C l] lines, (a) If we assume that the nar- 
row component is generated by noise and fit the 8727A line 
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Figure 12. The temperature sensitivity of the [C l] 9850A+ 
9824A/ 8727A ratio, in cases where the lines are generated from 
recombination, (s ohd lines labelled for Ca se A & B recombina- 
tion respectively, lEscalante &: Victoi]|l99oh and for cases where 
the lines are coUisionally excited (dot ted lines p l otted for various 
electron densities). Reproduced from llTu et al.l l(l99^. We have 
overplotted a box representing the likely margins of error of our 
[C l] ratio measurement and likely nebular temperatures. 



with one component, the [C l] ratio (9850 A+ 9824A/ 8727A) 
takes a value of around 8-9, which would usually indicate 
recombination as the origin of the lines, (b) Conversely, if 
we fit the 8727A line with a narrow and a wide component, 
we obtain a [C l] ratio of around 10-11, which is closer to 
the range of values which would indicate that the lines are 
coUisionally excited in nature. On Figure [TT] we have plotted 
a box representing the uncertainty range of the [C l] ratio 
and the likely te mperature range under consideration. 

Figure 2 of iLiu et al.l (j 19951 ). reproduced here as Fig- 
ure 1121 shows the sensitivity to electron temperature and 
density of coUisionally excited [C i] 9850 A 9824A/ 8727A 
(dotted lines). For Case B recombination (indicated by [C l] 
ratios of 8-9) the temperatu re dependence is very weak (solid 
lines) (|Escalante fc Victor|[l99Cli ) . It should be noted though 
that the 9824A and 9850A lines are emitted in a fixed in- 
tensity ratio (1:2.96), we observe these lines to have a ratio 
of 2.6 - suggesting that the uncertainty on the [C l] ratio is 
a little higher than that quoted. 

From Figure [12] we can see that at T=12,000 K the 
difference between the [C l] diagnostic ratios expected of 
Case B recombination and collisional excitation (using the 
lowest density track) is of the order 5-10. Using the quoted 
line ratio (2.96) we can find the [C l] ratio using either 9850A 
or 9824 A, the maximum [C l] ratio occuring when we use 
only the 9824 A line - yielding a value of around 11. The [C l] 
diagnostic ratio does not therefore unambiguously indicate 
the mechanism generating the lines. 

When the [C i] lines are observed as recombination lines 
they are usually thought to be emanating from relatively 
cold, dense, photo-diss ociation regions (PDR s) surround- 
ing the ionized regions. lEscalante et al.l (|l99ll ) showed that 
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Figure 13. Comparison of the velocity profiles of [O l] 6300A 
(dashed), [C l] 9850A (dot dash) and [O ll] 3726A (solid) in the 
UVES spectrum of NGC 3199. The narrow component of the [O l] 
6300A line arises from night sky emission. 



these lines can emanate from PDRs next to high density 
molecular material , which has be en detected in the vicin- 
ity of NGC 3199 ijMarstonllioOll ). This material could be 
the origin of the narrow component of the [C l] 8727A line 
(while providing a relatively small contribution to the ob- 
served 9824 and 9850A lines). The velocity offset of the peak 
of the narrow and broad components of [C l] 8727A is pos- 
sibly evidence of this. 

In Figure[T3]we show the velocity profiles of [O l] 6300A, 
[C i] 9850A and [O ii] 3726A. The [O i] and [O ii] lines 
are coUisionally excited in the nebula, with a characteristic 
temperature of 8000-10000 K. The line profiles are therefore 
thermally broadened (which we will use to derive tempera- 
tures in the next Section). The fact that the [C i] 9824 and 
9850A lines have very similar widths to those of the [O l] 
and [O ii] lines indicates that they are generated at similar 
temperatures. 



6.3 Line Broadening Temperatures 

The observed width of a spectral line depends upon several 
factors, primarily the nebular temperature, turbulence and 
expansion along with the instrumental line profile. Only the 
temperature broadening is a function of atomic mass, so if 
we can measure the widths of several lines from ions of dif- 
ferent atomic weight to sufficient accuracy, it is possible to 
derive the temperature of the line-emitting plasma. Histor- 
ically this has usually been performed using only a single 
pair of lines, [N ll] 6584 A and Ha, d ue to both their prox- 
imity and strengths (e.g. lDopitalll97 j ). The velocity profiles 
of [N ii] 6584A and Ha in our UVES observations are shown 
in Figure [m Given the resolution of the VLT/UVES obser- 
vations it is possible to accurately determine the widths of 
a large number of strong emission lines (including the tradi- 
tional pair) and to derive temperatures from more than one 
pair. 

The major obstacle which must be overcome to per- 
form this process is a detailed description of the intrin- 
sic line shapes. For coUisionally excited lines this is rela- 
tively simple: most have negligible intrinsic fine structure 



12 D.J. Stock et al. 

Table 9. NGC 3199 Nebular Diagnostics 



Ion 


Diagnostic 


Ratio 


Electron Density or Temperature 


[S II] 


6716A 


1.23 ± 0.01 


220 ± 20 cm~^ 


67314 






fO Til 


3729A 


1 22 + 01 


o^n + 20 cm~3 


37264 






[CI III] 


5717A 
57374 . 


744+0.069 
'^^-0.058 


500 ± 330 cm-3 


[S III] 


9069A+95314 


136 ± 2 


7600 ± 50 K 


63124 . 






[O III] 


5007A+4959A 


241 ± 7 


9600 ± 100 K 


4363A . 






[O II] 


3726A+3729A 


45.5 ± 0.5 


10350 ± 50 K 


7319A+7320A+7330A+7331A 






[O I] 


6300A+6364A 


> 51.5 


< 12250 K 


5577A . 






[S II] 


6717A+6731A 


13.9 ± 0.5 


9000 ± 1000 K 


4076A+40684 






[C I] 


9850A+9824A 


9.8 ± 0.6 




8 72 7 A 
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Figure 14. Observed velocity profiles of Ha (solid) and [N ll] 
6584A (dashed) in NGC 3199. The Ho line is broadened by 
fine structure components along with the thermal and turbulent 
broadening that also affect [N ll] . (Flux coordinate has been nor- 
malised to allow better comparison.) 



Table 10. Temperature determinations from line broadening in 
NGC 3199 



Species 


Line (A) 


FWHM (k 


m 






± 


Hi 


6562.80 


27.113 


± 





045 






[N II] 


6548.10 


13.411 


± 





055 


11580 


72 


[N II] 


6583.50 


13.505 


± 





050 


11520 


66 


[O II] 


3726.03 


13.250 


± 





062 


11570 


81 


[O II] 


3728.82 


13.373 







062 


11490 


80 


[S II] 


6716.44 


12.629 







087 


11560 


120 


[S II] 


6730.82 


12.597 


± 





080 


11570 


110 


[O III] 


4958.91 


15.975 


± 





054 


9730 


52 


[O III] 


5006.84 


15.840 







047 


9820 


47 


[S III] 


9068.90 


14.250 







142 


10570 


151 


[S III] 


9531.00 


14.695 


± 





175 


10290 


175 


[Ne III] 


3869.06 


15.036 


± 





086 


10260 


86 


[Ne III] 


3967.79 


14.578 







234 


10570 


241 


[Ar III] 


7135.90 


14.267 


± 





075 


10510 


82 


He I 


6678.15 


18.739 


± 





219 


9380 


157 


[Ci] 


9850.26 


12.152 


± 





230 


12490 


330 


[O I] 


6300.30 


12.175 







345 


12480 


510 



broadening and all broadening in excess of the instrumental 
broadening must be a combination of thermal and turbulent 
broadening. The situation is rather different for recombina- 
tion lines, which are often comprised of several very closely 
spaced fine-structure components. The first work done to de- 
scribe the fine structure components of the Ha line in this 
context found that the "natural" profile of the line could 
have a significant effect on t he derived electron tempera- 
ture (|Dvson fc MeaburnllToTll ). In our analysis we shall use 
the fin e structure correction term provided bv lClegg et al.l 
l|l999l) [Equation 9] for Ha, which follows the same form as 
those derived from earlier formulations while introducing a 
correction factor 5 which combines with the Ha FWHM in 
quadrature : 

FWHM(Ha)^ - 5'^ - FWHM([N ii])^ 

8fc[^^ ^] X ln(2) ^ ' 

where Te is the electron temperature, FWHM (Ha) and 
FWHM(N ll) are the FWHMs of Ha and [N ll], tuh and ruN 
are the atomic weights of hydrogen and nitrogen respectively 



and k is the Boltzmann constant. An explicit equation for 
calculating 5 is not quoted, b ut a table of values suitable for 
interpolation was included bv lClegg et al] l| 19991 ). 

We can generalise Equation [T] to apply it to other 
species, replacing FWHM([N ll]) and mjv with values ap- 
propriate for the line and species we wish to use. 

FWHM(Ha)^ - 5'^ - FWHM([X])^ 

^] X ln{2) ^ ' 

The results of applying Equation[2]for the strongest col- 
lisionally excited lines with Ha appear in Table [Tol Singlet 
He I lines can be included as they have no fine structure 
components. The only singlet He I line found with sufficient 
strength to measure its width accurately was He I 6678.2 A. 
We have not included temperatures derived from [C l] or 
[O l] in the plots as they are potentially unreliable. Equa- 
tion[2]is only valid if we believe the lines are being generated 
by the same gas. In the case of [C l] and [O l], this may not 
be the case as they are likely to be generated in the outer 
regions of the nebula. 

In Figure 1151 we show the line broadening tempera- 
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Table 11. Sources of Atomic data 



ID 



20 25 30 35 

Ionization Potential. eV 



40 



45 



Figure 15. NGC 3199 line broadening temperatures for various 
ions relative to Hcj, versus the ionization potential to obtain the 
parent ion. 
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Figure 16. Comparison of line broadening temperatures T(LB) 
with forbidden line ratio temperatures T(FLR). The blaclc solid 
line corresponds to T(LB) = T(FLR). 



ture versus tlie ionization potential for each species. The 
species with higher ionisation potentials generally have lower 
temperatures by at least 1000 K. This trend could be at- 
tributable to an effect known as "radiation hardening" (see 
lTielens|[2005l . Section 7.3.2). The ionization cross section of 
H I is proportional to , meaning that the lowest energy 
photons are absorbed by the inner regions of the nebula. 
Hence as radius increases the average energy and heating 
effect of the remaining ionizing photons also increases. For 
a radiation bounded nebula, these outer regions will also be 
where the degree of ionization decreases. 

Figure [16] shows the comparison between temperatures 
calculated using the line broadening method with those 
determined using diagnostic ratios. In most cases the line 
broadening temperature is higher than the diagnostic ratio 
temperature, [O ill] showing the closest agreement. 

6.4 Nebular Abundances 

For each nebula, diagnostic uncertainties were computed via 
a novel Monte-Carlo technique. The "NEAT" (Nebular Em- 
pirical Analysis Tool) software, developed at UCL, allows 
an exploration of the parameter space defined by the uncer- 
tainties associated with nebular line measurements. A full 
description of this code and its capabilities will be forth- 



lon 



Reference 



0° 

C1++ 

QO 

0+ 

S++ 

All 

others 



jPequignot fc Aldrovandil lll976l). Johnson et al.l 
tip7), Nussbaumer fc Ruscal lll979F 



iMendoza fc ZeippenI lll982bl ). iButler fc ZeippenI 

ii9gg|), 

Berring ton fc Burkj lll98ll) . iBerringtonI l ll98Sf) , 
Baluia fc Zeippe nl lll988l') 
PradhanI lll976l'). IZeippenI lll982l) 



PradhanI II1976I). IZeipnenI II1982I) 

IVIendoza fc ZeippenI il982a^ , llVlendozal lll983h 



Landi et al.l ||2006| ') 



coming (Wesson et al., 2011, in prep). Briefly, NEAT reads 
in a list of measured spectral lines (their wavelengths, rel- 
ative intensities and uncertainties), and generates a "possi- 
ble" line list by randomly offsetting the line intensities using 
the quoted uncertainties as limits. It then performs a full 
empirical analysis on this "possible" line list, calculating all 
diagnostic ratios and plasma parameters, along with a full 
abundance analysis using a three-zone ionization model by 
default. It repeats this process as many times as the user 
requires, building up a picture of the uncertainties associ- 
ated with each diagnostic, abundance, etc, as the number of 
data points increases. One then bins these values to obtain 
an uncertainty distribution. NEAT was used to generate all 
of the values and uncertainties in Table [5J along with the 
abundance derivations described later. 

The NEAT code uses the same statistical equilibrium 
methods as the earlier EQUIB code (developed at UCL by 
S. Adams and I. Howarth) to derive abundances for CELs, 
solving the equation: 

N{ion) _ ^ I(line) 



X A : 



(3) 



N{H+) I{HI3) ' 

where A is related to the emissivities of the two ions 
via the following relationship: 



e(line) n(ion)nee 
A oc -f oc 



(4) 



e{Hp) " N{H+)neT, 

The A value for the specific ion at given Ue and Te is cal- 
culated and then used to derive j^^^^"^^ , the linear abundance 
of the specific ion as indicated by the particular line chosen. 
NEAT derives uncertainties on the abundance determina- 
tions in the same manner as described for the temperature 
and density diagnostics. The provenance of the atomic data 
used is shown in Table [TT] 

In the simplest cases, only one of each type of diagnos- 
tic was available per nebula e.g. BAT99-38, where Te([S ill]) 
and ne([S ll]) were available. This led to a very simple abun- 
dance determination where the same temperature and den- 
sity was used for every ion. This was the situation for each 
of the three LMC WR nebulae. 

For NGC 3199, multiple density and temperature diag- 
nostics were available. Each density diagnostic was found to 
place the nebula at the diagnostic's low density limit. The 
forbidden line ratio temperature diagnostics indicated that 
a multi-zone ionisation model may be appropriate, in that 
lower ionisation species (e.g. [O ll] and [S ll]) gave higher 
temperatures than those from moderate ionisation species 
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(e.g. [O III] and [S in]). This trend was corroborated by the 
line broadening temperature measurements (see Figure [15)) . 
Ultimately it was decided to use a two zone model where 
we used the average of the [O ll] and [S ll] forbidden line 
ratio temperatures for the low ionisation species, the [S ill] 
temperature for [S ill], and the [O ill] temperature for all 
other ions. 

The results of the abundance calculations for the LMC 
WR nebulae and NGC 3199 are listed in TablefT2lusing ICFs 

1 Ionisation Correction Factors) from iKingsburgh &: Barlow! 
1991). In Table [T3l we list the N/0 abundance ratio in each 
of our nebulae along with some comparison objects. 

The remaining three Galactic nebulae in our 
NTT/EF0SC2 sample have no observed temperature 
diagnostic. As discussed earlier in this Section, the abun- 
dances from CELs depend on a Boltzmann factor and as 
such are sensitive to the adopted electron temperature. 

Previous authors have adopted temperatures for some 
of these objects based on matching the oxygen abundance 
to those of H II regions in the same galaxy. However we 
do not think this approach is appropriate in this case, as 
the oxygen abundance in WR nebulae could easily be de- 
pleted relative to Galactic H ll regions due to CNO cycling 
by the massive star prior to it becoming a WR star. One so- 
lution might be to assume a nebular temperature that yields 
Galactic H ll region neon abundances, as neon should not 
have been astrated by stellar processing. Unfortunately the 
higher reddening in the Galactic plane meant that we had no 
detections of [Ne ill] 3868, 3967A for ROW 58 or the nebulae 
around WR 8 and WR 16, and hence no neon abundance 
for these objects. 

However we can infer that the abundance pattern in 
these nebulae is not the same as in H ll regions, or indeed in 
the other WR nebulae in our sample, by simple inspection of 
the detected lines. WR 16 and ROW 58 in particular display 
[N ll] lines comparable in intensity to Hq, indicating likely 
nitrogen enhancement. 

In the absence of a temperature diagnostic for these 
nebulae, we employed a different approach to finding one 
important diagnostic - log(N/0). We know, from our ob- 
servations and from previous authors, that WR nebulae are 
likely to possess electron temperatures in the range 7000- 
12000 K. 

The EQUIB code, as mentioned previously, can be used 
to generate the appropriate line emissivities e across this 
range (using the atomic data listed in Table [TT|). We can 
then plot N"'"/0''' as a function of electron temperature Te 
as: 



I([N ii]6584A) 
I([0 ii]3727 + 3729A) 



X E 



where: 



E = 



[O II13727A 



+ e ° (To 

[O II13729A^ 



e °(Tc 
[N II]6584A^ 



(5) 



(6) 



Where I(line) is the line intensity relative to H/3 = 100 
and e(Te) is the line emissivity at temperature Te. This 
/O'^ relationship is shown in Figure [T71 for ROW 58 and 
the nebula around WR 8. As expected, the ratio is not a 
strong function of temperature and therefore provides a rea- 
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Figure 17. Range of possible log.^^^ values over the range 7000- 
12000 K for the nebulae around WR 8 and WR 40 (RCW 58). 



Table 13. N/O ratios 



Nebula 



log(N/0) 



BAT99-2 
BAT99-11 
BAT99-38 
LMC H II region 



^+0.07 
-0.10 
7+0.08 
-0.12 

-1.48 ± 0.05 



-1.52: 
-1.67^ 



-1.49 



WR 8 
WR 16 
NGC 3199 
RCW 58 

Galactic H II region 
Solar" 



> -0.4 

> -0.1 
-0.90 ± 0.01 

> 0.1 
-0.91 
-0.86 



'': lAsDlund et al.l l|2009h 



sonable constraint upon the value of log^. 
degree of nitrogen enrichment within each nebula. 



and hence the 



As mentioned earlier, we used the icf scheme of 



IKingsburgh fc Barlow! l|l994 ). which specifies identical icfs 
for N"^ and O""". This means that the N"^/0^ ratio that we 
calculate should be the same as the overall N/O ratio. The 
range of N/O ratios derived for RCW 58 and for the nebula 
around WR 8 will be discussed in Section [7] 



7 NOTES ON INDIVIDUAL OBJECTS 
BAT99-2 

The EF0SC2 optical spectrum of the nebula around BAT99- 
2 (WN2) i s in good agreement with that presented by 
iNaze et al.1 (|2003l l. The nebula is highly ionized which man- 
ifests itself via the presence of [Ar iv] lines along with He ll 
4686 A, which was only detected in BAT99-2. 

The nebula around BAT99-2 does not display any sig- 
nificant under- or over- abundances relative to LMC H ll 
regions (see Table [12]) . In g eneral the abunda nces we derive 
agree with those derived bv lNaze et aD (|2003l l and we agree 
with their conclusion that BAT99-2 does not contain any 
abundance anomalies. 
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Table 12. Total Abundances 



BAT99-2 



BAT99-11 



BAT99-38 LMC Hll Region" NGC 3199 Milky Way Hll region"" 



He+ / H+ 

He++ / H+ 

12 + log He / H 



0.042 ± 0.007 
0.066 ± 0.002 
11.04 ± 0.02 



0.102 ± 0.015 0.083 ± 0.014 



11 Ql+o-o'' 



10.92«;?« 



10.98 



0.088 ± 0.006 
10.95 ± 0.01 



10.98 



CO / H+ (xlO'^) 
12 + log C°/H+ 



5.97 ± 0.20 
5.78 ± 0.20 



N+ / H+ (xlO^) 
N ICF' 

12 + log N / H 



1^+0.50 
-0.30 
,+0.07 
-0.11 



4.65: 
6.61"* 



Q rjo-|-0.23 

^•^"-0.28 
1 cQ+0.25 
^■"^'^-0.07 

6.56 ± 0.06 



2.25 ± 0.20 

--I-0.50 
-0.05 
1-1-0.14 
-0.06 



1.45: 

6.54"* 



2.31 
6.92 



21.8 ± 0.6 
2.67 ± 0.02 
7.77 ± 0.01 



3.87 
7.49 



0° / H+ (xlOS) 
0+ / H+ (xlO^) 
0++ / H+ (xlO^) 
O ICF* 

12 -f log O / H 
12 -I- log 0°/H+ 



-1-0.04 
0.05 
-1-0.40 
0.30 
-1-0.17 
0.08 

12 ± 0.03 



2.80 
4.10 
1.83 



10 JO+O'^'O 
lU. (U_Q 20 
1+2.7 
-0.07 



5.90^ 



1 



8.23 



-0.06 
-0.02 



^■^^-0.37 
1+4.00 
-0.85 
1 

+0.14 
0.05 



3.20"* 



1.01 



8.41 



1.86 ± 0.09 
17.3 ± 0.8 
28.9 ± 0.8 
1 

8.66 ± 0.02 
7.28 ± 0.02 



1 

S.55 



log N / O 



-1.52 



+0.07 
0.10 



-1.67 



+0.08 
0.12 



-1.48 ± 0.05 



-1.49 



-0.90 ± 0.01 



-0.91 



Ke++ I H+ (xlO-"^) 1.13 
Ne ICF'' 3.07 
12 -I- log Ne / H 7.55 



+0.13 
0.10 
+0.33 
0.20 
+0.04 
0.03 



0.20^ 



1+0.12 
-0.03 

2.60 ± 0.88 

3+0.10 



1.83" 



•0.11 



0.44"* 



1+1.10 
-0.25 

2.60 ± 0.96 



7.39 



+0.24 
-0.40 



1.76 
7.62 



10.1 ± 1.0 
1.60 ± 0.01 
8.23 ± 0.02 



1.39 
7.99 



S+ / H+ (xlO") 
S+"^ / H"^ (xlO^*) 
S ICF' 

12 -I- log S / H 



0.46"* 



0.50 ± 0.04 

,+0.04 
-0.03 

1.25 ± 0.02 
6.08 ± 0.02 



0.45 
1.07 
1.02 
6.20 



0.07 
0.06 
+0.25 
0.09 
0.01 
0.01 
+0.06 
0.03 



1+0.08 
-0.06 
1+0.50 
-0.25 

1.01"^°°2 



0.59: 
1.02"* 



6.21 



0.01 
+0.13 
0.07 



6.70 



2.10 ± 0.06 
15.6 ± 0.50 
1.10 ± 0.01 
7.29 ± 0.01 



1.19 
7.00 



Ar++ / H+ (xl06) 
Ar3+ / H+ (xlO^*) 
Ar4+ / H+ (xlO^*) 
Ar ICF' 
12 -I- log Ar / H 



0.22"* 



1+0.03 
-0.01 

0.68 ± 0.01 
0.14 ± 0.04 

1 



6.02 



+0.04 



0.62 



-0.20 
-0.09 



1.87 
1.08 



+0.11 
0.07 



0.51 



+0.42 
0.11 



1.87 



6.06 



+0.21 



1.79 

6.30 



2.37 ± 0.05 



1.87 
6.65 ± 0.01 



1.35 
6.31 



a: iTsamis et all l|2003f l (LMC Nl lB), (M 17) 

h: Using the lonisation Correction Factor scheme of lKingsburgh fc Barlowl ((iSSJ) 



BAT99-11 

The nebula encircling BAT99-11 (WC4) displays a lower 
degree of ionization which is more typical of WR nebu- 
lae despite having a relatively high excitation (WC4) cen- 
tral star. The [Ne ill] and [Ar iv] lines seen in the neb- 
ula around BAT99-2 are not observed. Instead we see pro- 
nounced higher order Balmer series lines along with several 
He I lines, most prominently He I 447lA. 

The nebular abundances display no signs of stellar pro- 
cessing relative to H ll region LMC NllB (see Table [T2)) . 
Since BAT99-11 is a WC4 star it had been anticipated that 
the nebula might contain an ejecta component. Based on the 
lack of abundance anomalies, the nebula around BAT99-11 
is a wind-blown shell. 



BAT99-38 

The spectrum of the nebula near BAT99-38 (WC4-)-0) is 
very similar to that of the nebula around BAT99-11, albeit 
roughly a factor of five fainter in surface brightness. Many 



He I lines are present, including 7065A, 6678A, 447lA and 
several higher order Balmer lines are present blue-wards of 
4000A. 

The nebula near BAT99-38 displays no abundance en- 
hancements relative to H ll region LMC Nil (see Table [12)) . 
Based on our results and the fact that the star is superposed 
on the arc of nebulosity (Figure[3ll, it is dubious whether the 
nebula near BAT99-38 is actually related to BAT99-38 at all. 



WR 8 

The nebulosity surrounding WR 8 (WC7/WN4) presents 
clear evidence of nitrogen enrichment in the form of strong 
[N ii] 6548, 6584A emission lines relative to Ha. There is also 
some evidence of WR features scattered into the slit, since 
a wide bump at around 4650A corresponds to a known WR 
emission feature. 

From Figure [T71 we can see that the range of log(N/0) 
indicated, -0.4 < log(N/0) < 0.1, is well above the N/0 ra- 
tio found for the M17 Galactic H ll region (-0.91). Even the 
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lower bound, -0.4, is a factor of three above the Galactic ISM 
value, indicating a degree of nitrogen enha ncement relative 
to oxy gen. This supports the conclusion of lStock fc Barlowl 
1I2OIOI ) that this nebula has a significant stellar ejecta com- 
ponent. 



ICF. However, if we assume that NGC 3199 has a simi- 
lar to tal log(C/H) value to the Gala ctic H 11 region Orion 
(8.37; ISimon-Diaz fc Stasiiiskal I2OIII '). it would imply that 
only 0.25% of the total gas-phase carbon was in the neutral 
state. 



WR 16 



RCW 58 (WR 40) 



The nebulosity around WR 16 (WN8) displays a sparse, 
but interes ting spectrum. We no ted earlier that previous ob- 
servations (|Marston et al.iri999l ) had not detected the [S 11] 
doublet, a result repeated in our spectrum. The spectrum 
also contains the strongest [N 11] lines in our sample along 
with no trace of [O ill], [O 11] or [O l] lines, making a direct 
N/0 determination impossible but implying a high N/O ra- 
tio. This kind of composition is precisely what one might 
expect from stellar outflows of a WN type WR star. There 
is a strong sign of helium overabundance, with He I 6678A 
present and stronger than in all but the nebula around WR 
8. Given the non-detection of the [O 11] doublet at 3727A 
we cannot employ our 'bootstrapping" approach to calcu- 
late an approximate N/O value for the nebulosity, except to 
say that it must be strongly nitrogen enriched since if we 
assume a conservative upper limit for the [O 11] 3727A in- 
tensity, for example 'i^^^^^™^ ^ q 5 (as is the case for WR 
8 and RCW 58), then we find -0.1 < log(N/0) < 0.3 for the 
temperature range 7000-12000 K, a factor of 10 higher than 

the solar N/O valu e. 

iMarston et ai] (|l999l ) suggested that the lack of de- 
tectable [S 11] and [O ill] lines in the spectrum of the WR 16 
nebula was evidence of an extremely high nebular density 
(> 10000 cm"^). The critical densities of [S 11] and [O ill] do 
indeed lie in this range, however we attribute the weakness 
of the [O III] lines to the low effective temperature of the 
WN8 central star. 



NGC 3199 (WR 18) 

NGC 3199 displays very similar abundances to the Galactic 
H 11 region M17 (see Table [12}. We conclude that the nebu- 
losity is mainly a swept-up ISM shell. Previous observations 
of this nebu la have der i ved an expansi on velocity of around 
15 km s"^ (|Chul ll982': 'M arstonll200ll ). This value is simi- 
lar to the velocity separation of the main nebular emission 
components in our UVES spectrum (see Figure [TO)). How- 
ever, the low intensity of the red-shifted component suggests 
that it is unlikely to have the same intrinsic luminosity as 
the blue-shifted component, as the mass of internal dust re- 
quired to produce the required extinction (A(Ha) > 4mags) 
is unphysical. By inspection of Figure |S] it is obvious that 
NGC 3199 is not spherically symmetric, its shape suggesting 
a swept up origin, in which case the relative intensities and 
velocities of the emission components most likely are due to 
variations in the local ISM density. 

Following the discussion in Section 16.11 the dominant 
excitation mechanism for the observed [C l] lines in NGC 
3199 appears to be coUisional excitation of C° in the 
outer regions of the nebula. The C''/0'' abundance de- 
rived, log(C°/0°) = -1.49, would require detailed photoion- 
ization modelling to convert into log(C/0) as the physical 
and ionization structure of the nebula precludes a simple 



RCW 58 around the WN8 star WR 40 displays a spectrum 
with much in common with the other low-excitation Galac- 
tic WR nebulae in our sample. Strong [N 11] features are 
combined with weak [O ill] lines and a clear He I spectrum. 

From Figure [17] it is clear that RCW 58 has a very high 
measured log(N/0), lying in the range 0.1 < log(N/0) < 
0.55. The lower limit on this value (assuming a factor of two 
uncertainty) is still very high. This implies that the material 
possesses a CNO-cycle processed stellar ejecta component. 

This fi nding agrees wi th pr evious results, h o wever the 
analyses of iKwitterl l| 19841 ') and iRosa fc Mathij l|l990h re- 
lied upon the method of tuning the electron temperature 
until the "correct" oxygen abundance was achieved (in this 
case the Galactic ISM oxygen abundance). This methodol- 
ogy may be inappropriate for nebulae which are believed 
to contain processed material, as the true oxygen abun- 
dance is likely to differ from the ISM oxygen abundance. 
iKwitteJ ||1984 ) found log(N/ 0) = -0.42 usirig an e lectron 
temperature of 7500 K while iRosa fc Mathid (Il990h found 
log(N/0) = -0.3, however lRosa fc Mathid (|l990l ') did not in- 
clude a list of line detections, or finding charts indicating 
which part of RCW 58 they observed so it is impossible to 
draw further comparisons. In both cases the difference in 
derived log(N/0) stems from the difference in the strength 
of the [O 11] 3727, 3729 A doublet. We found a dereddened 
[O 11] doublet intens ity of ~ 42 (relative to H/3 = 100), 
whil e iKwitterl lll984h quote 172, a factor of four greater. 
The iKwitter l|l984h observations were not sky-subtracted, 
which may go some way to explain this discrepancy. 



8 CONCLUSIONS 

We found evidence of significant nitrogen abundance en- 
hancements in the nebula around WR 8, a newly discov- 
ered nebula which wc had morphologically categorised as 
an ejecta type nebula (S tock fc Barlow. ,201Q ). The Galac- 
tic nebulae, with the exception of NGC 3199, were found to 
have significantly enhanced N/O ratios compared to M17. 
The nebulosity around WR 40 (RCW 58) was the most en- 
riched (log(N/0) > 0.1), followed by that around WR 16 
(log(N/0) > -0.1), both of which have WN8h exciting stars. 
The nebula around WR 8 (WN7/WC4) is also significantly 
enriched in nitrogen (log(N/0) > -0.4). 

The N/O ratios derived for all of the LMC nebulae were 
consistent with LMC H 11 region values (log(N/0) = -1.49). 

Our very high resolution UVES data allowed us to de- 
termine the NGC 3199 spectral line widths very precisely. 
In combinat i on wi th correction factors for Ha published by 
IClegg et al.1 (|l999l ). these widths were used to measure the 
nebular temperature by determining the degree of thermal 
broadening thought to be affecting each line. The tempera- 
tures determined using this method for different ions were 
systematically lower than those determined using empirical 
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line ratio temperature indicators for the same ions. Higher 
electron temperatures were found for low ionization ions - 
an effect we attribute to radiation hardening. 

A secondary goal of this project was to determine 
whether carbon abundances could be derived for nebulae 
around carbon rich WR stars. Insufficient spectral resolu- 
tion precluded us from detecting [C i] lines from the WC 
nebulae which we targeted. However we have shown that 
the far-red [C i] lines could be detected from NGC 3199 
in a relatively modest amount of time with a suitably high 
resolution instrument. 

We did not detect any heavy element recombination 
lines in the UVES spectra of NGC 3199. Since these hncs lie 
mainly in the blue this was likely due to the high reddening 
towards this Galactic plane target. A similar set of UVES 
observations of suitable nebulae in the LMC could poten- 
tially yield detections of the far-red [C i] lines and the C ii 
4267A recombination line, which could allow a quantifica- 
tion of the carbon abundances in WR nebulae. 

The lack of appropriate temperature diagnostics for 
several nebulae rendered the derivation of elemental abun- 
dances for them impossible. The non-detection of the [O iii] 
4363A line in several Galactic WR nebulae was due to sev- 
eral effects: low spectral resolution, high reddening and low 
O^"*" abundances. Observations at higher spectral resolving 
powers, similar to those conferred by UVES, can enable the 
detection of crucial temperature diagnostic lines, along with 
the above-mentioned diagnostic lines of carbon. 
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